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We study the density-density correlation function of the dense random het- 
eropolymer solutions. We show that a phase transition is possible due to the het- 
erogeneity of polymers. We also show that the critical behavior of the system is 
described by the 0{N) model at A'' = 0. 



PACS number(s): 61.25.Hq 



I. INTRODUCTION 



Random heteropolymers are of great interests because of their wide occurence in nature, 
e.g. many biopolymers are random heteropolymers. It is well known that they have many 
unusual properties For example, when the temperature decreases they experience 

microphase separation and freezing phase transitions. The protein folding is also an example 
of such a property. However all of these properties refer mainly to a single heteropolymer. 
The properties of polymer solutions have been well studied [^,0. One of the well known 
properties is the screening of excluded volume interactions in polymer solutions. It leads 
to an exponential decrease of the density- density correlation function with distance. In two 
dimensions, some peculiarities appear in the behavior of the correlation function. It turns 
out that the excluded volume effects lead to a phase transition in a dense polymer solution 
0. The ordered phase is described by the low temperature phase of = vector model 
101 . One of the features of the ordered phase is that the correlation function increases with 
distance. 

Along with the single heteropolymers and polymer solutions, it is interesting to inves- 
tigate also heteropolymer solutions. In this paper this problem is investigated for the first 
time. We will see below that a phase transition is possible due to the heterogeneity of 
polymer chains. 

The paper is organized as follows. In Section II we derive an equation for the correlation 
function of a random heteropolymer solution. In Section III we develop a diagram technique 
to treat with the heterogeneity random field. A critical heterogeneity parameter is introduced 
in Section IV. In Section V we show that the crossed diagrams for self-energy are smaller 
than the non-crossed ones. Then we sum the non-crossed diagrams and improve the mean 
field approximation. In Section VI we present a field-theoretical approach and show that 
the problem is reduced to the 0{N) model at A^ = 0. Section VII concludes our results. In 
the Appendix, we justify the approximations done in Section II. 
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II. INITIAL RELATIONS 



The correlation function of a polymer solution is determined as follows 



where 



G{r) = -[< c(f)c(0) > -c' 
c 



(1) 



(2) 



is the microscopic concentration of monomers, Ri is the coordinate of the i — th monomer, 
< ... > means the thermodynamic averaging and c is an average concentration of monomers. 
Note that Eq.(0) implies a summation over all chains and monomers. Now let us derive an 
equation for the correlation function (^. It can be represented in the form [^] 



-G(f) =aJ Y[dRic{f)c{0)exp 



U_ 
T 



(3) 



where U is the conformational energy of polymers and A is a normalization constant 



IJ_ 
T 



^ = J Yl dRiexp 

The conformational energy consists of elastic and excluded volume parts 

^ = ^ 51 {^bn - Rbn-l) + 2 51 Vij^iRi " Rj)- 



(4) 



(5) 



Here a = 3,2 for three and two dimensions, respectively, the index b denotes a polymer 
chain and n is the number of a monomer in a chain, a is the average size of monomers, 
i = b,n and Vij are the random excluded volume constants. It is convenient in (Q) to go to 
integration on the monomer concentration c(r) instead of Ri. For this purpose one has to 
express the conformational energy U in terms of monomer concentration. It is well known 
that the elastic energy of a polymer solution, in the Gaussian approximation, can be 
represented in the form 

Uel BV 



T 



ctc 



(6) 



where B = l/12c, l/8c for three ^ and two dimensions, respectively, V is the volume of the 
system and is the Fourier transform of c(r) 

1 



eg = — / c{r)e'^^df. 



(7) 



Note that the approximation (j^) is justified in the concentrated solutions where the density 
fluctuations are small It is convenient to represent the random excluded volume constant 
in the form Vij = v + 6vij, where v is the average excluded volume constant and 6vij is the 
fluctuating part with the average < 6vij >= 0. Using Eqs.(|D and (|^) for the average part 
of excluded volume energy, one finds 

Jjav 



T 



^ y dr(?{f). 



(8) 



It is shown in the Appendix that if bvij are independent random variables with the variance 
< bvlj >= P] then the fluctuating part of excluded volume energy can be represented in 
the form 



where v\{f) is a Gaussian distributed random function with 5 correlations 



(9) 



w 



< vi{f)vi{r ) >= — 6{f — r), < vi >= 0. 



(10) 



Using Eqs.(P), (1) and @ for the conformational energy of a random heteropolymer solution, 
finally we obtain 

U U,i + u^i + Ul^ 1 



T 



T 



'^""^ ' -exc _ - ^ ^^^^^ _B^2y2 ^ ^ ^ ^^^^ 



(11) 



Now the correlation function ([1|) can be written in terms of a functional integral over c(r) 



cG{r) = A f Dc{r)c{r)c{0)exp 



f/[c(f)] 



(12) 



where 



Dcif) = n dc^,. 



(13) 



fc>0 



The symbol nfc>o means an integration over the independent components cg. Because the 
relation c| = c_g (this follows from the fact that c{r) is a real function) not all the com- 
ponents of Cj: are independent. Calculating the Gaussian integral over c{r) in Eq.(|12]), one 
finds 



cG{f)= -Ba'V' + v + vi{r) S{f). 



(14) 



It is easy to obtain the equation for G{'r) from Eq.(|T 

,vi(r 



G{r) = -—S{f). 
vc 



(15) 



Here ^ = (Ba^/vY^"^ is the correlation length of the polymer solution. 



III. IMPURITY DIAGRAM TECHNIQUE 

For bare correlation function with vi = 0, which corresponds to homopolymer solution, 
one has from (p!5|) 

^o(^) = / 2I, 2V (16) 

In the coordinate representation we have from Eq . (p^) 

,^ f dq expjiqr) 



Using the method of steepest descent in the integral of Eq. ([T7|) over the momentum 0, we 
can find the asymtotic behavior of correlation function on large distances 

It follows from Eq.(|18|) that the correlation function at large distances decreases exponen- 
tially with the correlation length ^. 
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Now we construct a perturbation theory with the heterogeneity parameter vi. It is 
convenient to deal with the dimensionless in ci = 2 correlation function G{r) = G{r)vc^'^ 
that satisfies the equation 



72 c~2 c^-2^l(^) 



V 



G{r) = -5{r). (19) 



In order to construct a perturbation theory in the heterogeneity parameter, we use the 
impurity diagram technique [|l^. It is easy to obtain the following expansion from Eg . (1191) 



G= = + E + Z Z + , (20) 



where the self-energy S is determined by the following irreducible diagrams 



(21) 



The heavy line in Eq.pO|) denotes the averaged correlation function, the thin one denotes 
the bare correlation function of Eq. (0) and the dashed one is the Fourier transform of 
the heterogeneity random field correlation function of Eq.(|l5|), C.^'^w'^ / {cv"^) . Summing the 
diagrams in Eq.(^), we have the following Dyson equation 



(22) 



It has a solution 
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where the self-energy is determined by diagrams of Eq. (|2T|). In the leading order, we get 
for the self-energy 

So(g) = ^ 7 / 7T%'^o(g - k) = ^ ^ ln(l + — ) for d = 2, and 

for d = 3. (24) 



1 TT 



The divergence of the integral in the upper limit in Eq.(^) is caused by the 6 feature 
of the excluded volume interaction of Eq.(|^). It is obvious that it should be smoothed on 
the scales less than v^^'^. Therefore we cut the upper limit of the integral in Eq.(p^ by this 
value. 

It follows from Eqs. (|23|) and (^^ , that the "mass" term in the denominator of correlation 
function of Eq. (123|) is renormalized as 

= - S. (25) 



It follows from Eqs. (p3D -(p5D, that the second term in Eq.( |2^) at some critical value of 
heterogeneity parameter w'^/v'^ can become equal to the first one. In this case the correlation 
length diverges and the exponential decay of the correlation function of Eq. (|18]) at large 
distances will be substituted by a power law. This is a second order phase transition caused 
by the heterogeneity of polymers. 

IV. CRITICAL HETEROGENEITY PARAMETER 

It is easy to obtain the critical value of the heterogeneity parameter g = w'^/v'^ from the 
condition of divergence of the renormalized correlation length. Using Eqs. (p^)-(|25|), one 
has 

= ln( l'+1v^) 

for d = 3. (26) 



1/^1/3 _7r/2^ 

Taking into account that almost always ^ ^ v^^^, we can simplify the expression of Eq. 
Dby 



vra^ 



4wln^ 



for d = 2, and 



where $ = ca^ is the fraction of monomers in the solution. It follows from Eg . (1271) that 
Qc does not depend on monomer concentration (only logarithmically in two dimensions). It 
is mainly determined by the stiffness parameter f/a^ {v/a'^ in two dimensions) of polymer 
chains. It is also obvious from Eqs.(^6])-(^) that the phase transition in heterogeneity 



parameter is more appropriately realized in a system of flexible chains v/a^ < 1 rather than 
in stiff v/a^ <^ 1 chains solution because the limit g ^ is easier reached in the first case. 
One can obtain from Eqs.(p3|)-(|25|) that near the critical point the renormalized correlation 
length diverges as 

U - a^^V (28) 
9c 

So, for the critical exponent of the correlation length in the phase g < gc we obtain the 



mean field result. This means that tree approximation that we used in Eqs. (|20|) , (|23D 



and (p^ for correlation function corresponds to the mean field theory. At the critical point 
9 = 9c one has power law behavior for correlation function on large distances instead of the 
exponential decreasing in g < gc case. In the next section we try to go beyond the mean 
field approximation. 

V. BEYOND THE MEAN FIELD APPROXIMATION 

Let us consider the second order diagrams in the diagrammatic expansion of the self- 
energy in Eq.(^T]). For the contribution of the second diagram in (^), one has 

^■M = / §p|G'(*. - «1G(tl + k\ - ,T (29) 

We take the correlation functions in ( p9D in the tree approximation 

G{k) = (30) 
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where the renormahzed correlation length is determined by Eq.(p5D and ^ oo at the crit- 
ical point g = Qc- Substituting Eq. ( pOD into Eq. (pQ]) and integrating over the momentums, 
we have 







Cf 2 ■ 


' 167r2 








' 87r3t;i/3 



V 



for d = 3. (31) 

When obtaining Eq.(^) we cut the upper limit of momentum integrals in Eq.(^) at 
and believe that ^ v^^'^ which is always correct in the critical region g g^.. Now 
consider the contribution of the third diagram in Eq.(^). It can be represented in the form 

S3(g) = C-^r J -j^G{h)G{h + h)G{h - q). (32) 

Substituting Eq. (pOD into Eq. (|3T|) and going to the dimensionless variables of integration, 
one has 

We are mainly interested in the large distance behavior of the correlation function. Our 
analysis of the integral in Eq.(Rl) at the limit q —>■ shows that 



-rW)2 a 



^3(0) ~ {^-^Yi^. for d = % and 



for rf = 3. (34) 

CV^ OTT'^ V^l-^ 

Comparing Eq.(p4Dwith the Eq.(|3TD one can see that in the critical region ^/j S> v^l'^ 
the crossed diagrams ,because the additional integration over angles, (see also ||10[) have an 
additional smallness on parameter v^f'^ /^r <^ 1. So, in the critical region the non-crossed 
diagrams are the dominant ones and we can sum them 



z= ■- ^ + ^ + '^—^ +... 
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The sum of the diagrams can be represented in the form 

So 



= T-TPT for d = 2, and 

= 1 4 2 ^ for d = 3. (36) 

The expressions are correct provided that ^~^w^^/j/47rcf ^ < 1. At the critical point — * oo. 
Therefore, even this improved approximation which corresponds the one-loop correction in 
field theoritical calculation still breaks down. To investigate the critical region g ^ Qc more 
carefully and to describe the ordered phase g > gc, we use a field theoretical approach in 
the next section. 



VI. FIELD THEORETICAL APPROACH 



It follows from Eqs.(|T9|)-(|25|) and (^) that the renormalized correlation function satisfies 
the equation 



- + m^ + h{r) G{r) = 5(f) 



(37) 



where = ^^"^{gc — g)/ gc is the renormalized "mass" term and h{r) = ^~'^vi{f)/v is the 
heterogeneity random field. To carry out the averaging over the realizations of random field, 
we use a functional integral to represent the correlation function: 

-i f Dip{f)(f{r)Lp{o) exp 5*0 



G(r) 



J Dip (r) exp So 



(38) 



where 



So 



dfipif) — + + h{f) ip{r) 



(39) 



Here (p{r) is a real field. Introducing the replica fields and using Eq. (0), one can represent 
the correlation function in the form 



G(f) 



N 



(40) 
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where 



S, 



N 



drifaif) — + + h{f) ipa{r)- 



(41) 



In Eq.(^), a summation over the repeated indexes is imphed. Carrying out a Gaussian 
integration over ipa in denominator of Eq. (PO]) , one has 



Y[ D(pa{r) exp Sn = det 



-27ri 



(42) 



Finally, taking into account Eqs.(0) and (^2|) we obtain the following representation for the 
correlation function 



G{r) = hm -z J Y[D(pair)(piir)(piiO) exp Sn- 



(43) 



Now using Eq. (^31) , one can carry out the Gaussian averaging over the random field h{r) 



G{r) = hm -z j [| Dv9„(f)v9i(f)v9i(0) exp S", 



(44) 



where 



S 



m 



/■ r 2 / ^l2 

^a{r} - ^ dr ip^{r} 



(45) 



Here H = ^^^w"^/ (cv^). So, we obtain a non-trivial field-theoretical model. Expanding exp 5" 
in heterogeneity parameter H one can reproduce the impurity diagrams in Eq. (plj). 

So we map the concentrated random heteropolymer solution problem onto the 0{N) 
vector model at = 0. This model is well studied and we can use it's results. Remind 
that the problem of self-avoiding walk of a single polymer chain is also mapped to this model 



11|. However there is an important difference between these two problems. In our case the 



ordered phase g > Qc (see below) also exists in contrary the former one where the criticality 
is reached in the limit of long chains L —>■ oo. 

Thus, we have the following behavior for the correlation function in the phase g < gc- 
When g <^ g correlation function at large distances exponentially decreases with some 
correlation length. At g ^ g^ correlation length diverges ^ ^ {gc — g)"^ ■ The value of 
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the exponent v is well known from the e — expansion and from the numerical calculations 
u ~ 3/5, 3/4, for three and two dimensions, respectively. At the critical point the correlation 
function has a power law behavior G{r) ~ rjj^g exponent rj is also well known 

for 0{N) model. This power law behavior is preserved also in the ordered phase in three 
dimenions. In two dimensions there are some peculiarities in the ordered phase. It follows 
from the preceding consideration that the phase transition is caused by the "massless" 
fluctuations. In the vicinity of the critical point we can simplify the action of Eq. f^SD by 
separating the "massless" fluctuations. For this reason we find saddle-point trajectories of 
the action of Eq. (^5]). Differentiating (^Sf), we obtain 



iV'^^p{r) - imV/3(r) - —^^{r) ^ i^l^{r) = 0. (46) 

^ 7 

We are looking for the homogeneous solution of Eq. (^61), f/sir) = const 

- im'^iff^ - —^p V'? = 0- (47) 



2 7 



This equation has two solutions 



y.^ = 0, /3,7 = l,....iV. (48) 

7 ^ 

It can be seen from Eq.(|48|) that the second type of saddle points lies in the complex plane. 
In the stationary phase method it is well known that in such cases one must transform the 
integration contour so that it passes through the saddle points (Fig. 1). We have in mind that 
the integration contour is transformed locally at each point r. Along the new contour if has 
the form = Lp^{l—i), where ip^ is real. Suppose that the fluctuating fields locally satisfy the 



extremum condition of Eq. (^S|). Making a change of variables ^a{^ = yf^'^/H(p^{l — i), we 
obtain the following expression for the effective action describing the "massless" fluctuations 

2 

Seff = -^j driV^.r, (49) 

where 

Ev^a(r) = l. (50) 

a 
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Here, for brevity, we omit the index R and neglect the constant terms in action that do not 
affect the critical behavior. So, the critical behavior of random heteropolymer solution is 
described by the non-linear a model at = 0. The coupling constant of the model of Eq. 
(^) is rn? /H = cC,'^{gc — g)/{ggc)- The region g < gc corresponds to the high temperature 
phase of = vector model |T2[. Note that our derivation of a model from the (v^^)^ model 



is correct in the high temperature phase g < gc- It is well known [P] that at the critical 
point the correlation functions of the non-linear cr-model and (v^^)^ theory are identical. 
However (v?^)^ model is renormalizable at d < A but cr model is renormalizable only at 
d < 2. So, in the critical domain for d = 2 one can use the non-linear a-model instead 
of the (y9^)^-model. The ordered phase g > gc will be described by the low-temperature 
phase of the non-linear cr-model at = 0. However the coupling constant in this case 
differs from that of Eq. (|50D . It will be proportional to the expectation value of the order 
parameter which is non-zero in the ordered phase 0. Note that this transition to the low- 
temperature phase of = vector model in two-dimension is caused by the heterogeneity 
of chains and not by dense polymer self-avoiding walks The low-temperature phase of 
A^ = vector model has many unusual properties. One of them is the increasing of spin-spin 
correlation function with distance In our case this means that in the phase g > gc the 
correlation function will increase with distance. The corresponding power index is 3/8 . 
In the dense polymer case such a behavior is explained by the repelling of the extremities 
of chains. To understand this result for the heteropolymer solution, consider the case when 
each chain consists of only two kinds of monomers (see Appendix). In this case there will 
be three excluded volume constants Vaa, Vab and Vbb. Suppose that the excluded volume 
constants between monomers largely differ from each other so that Vab, Vhb ^ Vaa (see also 
Appendix). In order to minimize the conformotional energy, type a monomer prefers contact 
with type a monomers. This means that the compensation of excluded volume interactions 
leading to the screening effect in polymer solutions will be violated in random heteropolymer 
solution. Moreover, the strong repulsion between different kinds of monomers will lead to 
the microphase separation like the analogous effect in single random heteropolymers (see 
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IH] and and references therein). Although these papers mainly dealt with the compact 
conformations of single random heteropolymers (which is the case f < 0, in our case always 
f > 0), nevertheless we think that these two phenomena are closely related. The microphase 
separation is the main reason of the rise of anti-correlations in heteropolymer solutions in 
two dimensions and power law decreasing of correlation function in three dimensions. 



VII. CONCLUSION 

We have considered the effect of random heterogeneity on the density-density correlation 
function of dense polymer solutions. It turns out that a phase transition is possible due to 
the heterogeneity of the polymer chains. The heterogeneity parameter is determined by the 
ratio of variance and average values of excluded volume constant. The critical behavior is 
described by the 0(N) vector model at = 0. The anti- correlation behavior of correlation 
function in two dimensions and the power law behavior in three dimensions in the ordered 
phase are associated with the microphase separation. 

This work was supported in part by the National Science Council of the Republic of 
China (Taiwan) under the grant no. NSC 89-2112-M001-005. 



Appendix: Random Excluded Volume Constant 



Using (^, one can represent the fluctuating part of excluded volume energy in the form 

^ = lJ drY.^v,Ar-Ri)Kr- R,)- (51) 

Suppose that 5vij are independent random variables with variance w. Then for the large 
number of monomers because of the central limit theorem, one has 



^ 5v,,5{r- R,) = v,ir)cir}. (52) 
j 

where Vi are random Gaussian variables with zero average and variance w. To prove Eq. 
(^), we calculate the averages and variances of left and right hand sides of Eq. (|52D . It is 
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evident that < Vi >= because < Svij >= 0. Now calculate the variance of the left hand 



side of Eq. (pi), 



< 6vij6{r- Rj)6vikS{f- R^) >=Y< Svij6vik > S{r- Ri)6{f- Rk). (53) 

jk jk 

Because 6vij are independent random variables, one has 

< SvijSvik >= Sjkw'^. (54) 
Substituting Eq.(p4D into Eq.(|^), we have 

J2 < Svijdvik > 6{f- Rj)6{f- Rk) = w'^c^ir). (55) 

jk 



It is evident that the variance of right hand side of ( p2D is also w'^c^{r). Substituting Eq . (^D 
into Eq.(^) for the fluctuating part of excluded volume energy, we have 

^ = ll dfc{f)Y^v,{R,)5{r - R^. (56) 

Let us prove that 

Y^ViS{f- Ri) = vi{r)c{r), (57) 

i 

where Vi{r) is a Gaussian distributed random function with 5-correlations. It is evident that 
< J2i Vi5{r — R) >= 0. Now let us calculate the correlator 

K{f, f^) = Y,< ViVj > 6{f- Rj)6{r - Rj). (58) 

ij 

The Vi are independent Gaussian variables, therefore < ViVj >= 6ijw'^. Substituting this 
into the Eq. (^8]) , we have 

K{f, r') =w^Yl ^^Ar- R:W - R,). (59) 

— * 

For the dense solution of heteropolymers one can go to the continious limit 5ij S{Ri — 

Rj)/c, 

,2 



in 

K{f,7') = — - RM^- RiW - Rj)- (60) 
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Using the properties of S functions, one has 

tip 

K{fy) = —5{r- r')c{r)c{f^). (61) 

So the function J^ViSlf— Ri) is equivalent to the function vi{f)c{r). Finally, for the random 
part of excluded volume energy, we have 

^'"'^ - ^ / dfvi{r)c\r) (62) 



T 2 
which coincides with Eq. (^. 

Now consider a special case when each chain contains only two types of monomers. 
Suppose that each heteropolymer consists of a and b type monomers with 1 — / and / 
fractions, respectively. Consider the case when / ^ 1 and Va <^ Vf,. Here Va = Vaa and 
fj, = Vah- As the effect of v^b will be proportional to <C 1 , we can ignore it. In such 
conditions the average excluded volume constant can be estimated as 

V={1- f)Va + M- (63) 

The variance (dispersity) of the excluded volume constant is estimated as 

W^ = f{l-f){Vb-Vaf. (64) 

The heterogeneity parameter for this type of randomness will be 

v' |(i-/K + Mf ^ 

Assume that fvi, <C (1 — f)va then we get from Eqs. (|63|) and (|65|) 

v^Va. f%. (66) 

Although fvb/va <^ 1, the heterogeneity parameter can be large 

9 • (67) 

Va Va 



because of Vb ^ Va- As can be seen from Eq.(|66D, in this particular case the average excluded 
volume constant is determined mainly by type a monomers and the dispersity is determined 
by type b monomers. 
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Figure Caption 



Fig. 1. Transformation of the integration contour in the complex ip plane. 



18 



FIGURES 

FIG. 1. 

A 




Fig.l Transformation of the integration contour in the 
complex 9 plane 
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